DRAFT VERSION AUGUST 27, 2008, ACCEPTED FOR PUBLICATION IN THE ASTROPHYSICAL JOURNAL LETTERS 
Preprint typeset using L^T^X style cmulatcapj v. 04/21/07 



FORMATION OF A QUASAR HOST GALAXY THROUGH A WET MERGER 
1.4 BILLION YEARS AFTER THE BIG BANG 

DOMINIK A. RlECHERS 1 ' 2 ' 5 , FABIAN WALTER 1 , CHRISTOPHER L. CARILLI 3 , FRANK BERTOLDI 4 , AND EMMANUEL MOMJIAN 3 
draft version August 27, 2008, accepted for publication in the Astrophysical Journal Letters 

ABSTRACT 

We present high-resolution Very Large Array imaging of the molecular gas in the host galaxy of 
the high redshift quasar BRI 1335-0417 (z = 4.41). Our CO(J=2^1) observations have a linear 
resolution of 0.15" (1.0 kpc) and resolve the molecular gas emission both spatially and in velocity. 
The molecular gas in BRI 1335-0417 is extended on scales of 5 kpc, and shows a complex structure. 
At least three distinct components encompassing about two thirds of the total molecular mass of 
9.2xlO lo M0 are identified in velocity space, which are embedded in a structure that harbors about 
one third of the total molecular mass in the system. The brightest CO(J=2^1) line emission region 
has a peak brightness temperature of 61±9K within 1 kpc diameter, which is comparable to the 
kinetic gas temperature as predicted from the CO line excitation. This is also comparable to the 
gas temperatures found in the central regions of nearby ultra-luminous infrared galaxies, which are 
however much more compact than 1 kpc. The spatial and velocity structure of the molecular reservoir 
in BRI 1335-0417 is inconsistent with a simple gravitationally bound disk, but resembles a merging 
system. Our observations are consistent with a major, gas-rich ('wet') merger that both feeds an 
accreting supermassive black hole (causing the bright quasar activity), and fuels a massive star burst 
that builds up the stellar bulge in this galaxy. Our study of this z>4 quasar host galaxy may thus be 
the most direct observational evidence that 'wet' mergers at high redshift are related to AGN activity. 
Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Great progress has been made in recent years both ob- 
servationally and theoretically to further our understand- 
ing of galaxy formation and evolution from the early 
epochs of galaxy formation to the present-day universe. 
One basic prediction of cosmological simulations is that 
during the early epoch of hierarchical galaxy formation, 
some of the most massive galaxies are al ready formed 
in major merger events (e.g., Springel etal. 120051 ). These 
mergers are believed to commonly trigger both AGN and 
starburst activity in such early systems, w hich i s regu- 
lated via AGN feedback (e.g., Hopkins et al. 120051k Such 
feedback may be responsible for the present-day "Mbh~ 
a v " relation between black hole mass a nd bulge veloc- 
ity di spersion (Ferrarese & Merritt [2000l Gebhardt etal. 
l2000l) . 

Studies of molecular gas (the requisite material to fuel 
star formation) in young galaxies are an essential ingre- 
dient to understand their physical properties in more de- 
tail. Molecular gas (typically CO) has been detected 
in ^50 galaxies at z>l to date, revealing large molecu- 
lar reservoirs of >10 10 M Q in most cases (see Solomon 
& Vanden Bout 120051 for a review). However, these 
studies rely almost exclusively on the integrated prop- 
erties of the line emission, as the molecular reservoirs 
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in these distant galaxies are difficult to resolve. To 
date, only the z=4.69 and z=6.42 quasars BR 1202-0725 
and SDSS J1148+5251, two of the most distant gas-rich 
galaxies could be resolved in molecular gas emission with- 
out t he aid of gra vitati onal magnific ation (Omont et al. 
1996; Carilli etal-HH Walter etal.[200l). 

In this letter, we report on high angular resolution 
(0.15"; 1.0 kpc) Very Large Array (VLA) 6 observations 
of molecular gas in the host galaxy of BRI 1335-0417, a 
dust-rich, optically identified quasar at a redshift of 4.41, 
corresponding to only 1.4 Gyr after the Big Bang. Op- 
tical imaging with the Hubble Space Telescope reveals 
a single point source without any evid ence f or gravita- 
tional lensing (Storrie-Lombardi et al. IT996T ). We use 
a concordance, flat ACDM cosmology throughout, with 
ff n = 71 km s'^pc' 1 , ft M =0.27, and fl A =0.73 (Spergel 
etal. I2007T) . 

2. OBSERVATIONS 

We observed the CO(J=2^1) transition (^ rcst = 
230.53799 GHz) towards BRI 1335-0417 using the VLA 
in B configuration on 2005 April 11, 13, and 22, and in 
C configuration on 2005 September 16. The total on-sky 
integration time amounts to 29 hr. At z = 4.4074, the 
line is redshifted to 42.6338 GHz (7.03mm). Observa- 
tions were perfo rmed in fast-switching mode (e.g., Car- 
illi & Holdaway |l999) using the nearby source 13515- 
01513 for secondary amplitude and phase calibration. 
Observations were carried out under very good weather 
conditions with 25 antennas. The phase stability in all 

6 The Very Large Array is a facility of the National Radio As- 
tronomy Observatory, operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 
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Fig. 1. — VLA map of the CO(J=2^1) emission to- 
ward BRI1335-0417 (integrated over the central 43.75 MHz, or 
308kms — 1 ), imaged using the combined B and C array datasct 
and natural weighting. Contours are shown at (-3, -2, 2, 3, 4, 5, 
6, 7)xcr (ler = 50/iJy beam -1 ). The beam size (0.23 "x0.18 ") is 
shown in the bottom left corner. 

runs was excellent (typically <15° for the longest base- 
lines). The phase coherence was checked by imaging 
the calibrator source 13569+02144 with the same cal- 
ibration cycle as that used for the target source. For 
primary flux calibration, 3C 286 was observed during 
each run. Given the restrictions of the VLA corre- 
lator, two 25 MHz intermediate frequencies (IFs) with 
seven 3.125 MHz channels each were observed simultane- 
ously centered at the CO( J=2— >1) line frequency, leading 
to an effective bandwidth of 43.75 MHz (corresponding 
to 308 kms -1 at 42.6 GHz). This encompasses a ma- 
jor fraction of the CO line width as measured in the 
CO(J=5-> 4) tra nsition (420±60 kms -1 FWHM, Guil- 
loteau et al. 119971) . but does not cover the line wings and 
the continuum. Earlier observations set a 2 a li mit o f 
240/xJy on the continuum emission (Carilli etal. Il999f ). 
in agreement with an expected flux of ~34 /zJy as derived 
from the continuum spectral energy distribution of this 
source. 

For data reduction and analysis, the AIPS package 
was used. Two data sets were created for the final 
analysis. The first dataset includes both the B- and 
C-array data, and is imaged using natural weighting. 
A velocity-integrated CO( J=2—*l) map of this dataset 
is shown in Fig. [T] The synthesized clean beam has 
a size of 0.23"x0.18" (1.6x1. 2kpc). The final rms 
over the full bandwidth of 43.75 MHz (308 kms -1 ) is 
50 /iJy beam -1 . To boost the resolution, the second 
dataset includes B array data only, and are imaged 
using robust (i.e., intermediate between natural and 
uniform) weighting (Fig. [J), achieving a resolution of 
0.16"x0.14" (l.lx0.95kpc). This results in an rms 
of 58 /xJy beam -1 per channel. In Fig. [3l seven veloc- 
ity channel maps (6.25 MHz, or 44 kms -1 each) of the 
CO( J=2— >1) line based on the combined BC-array data 
are shown. A Gaussian taper falling to 30% level at a 
baseline length of 800 kA (longest baseline in the B array 
data: ~1.6MA) was applied to the UV data, resulting in 
a resolution of 0.32 "x 0.30" (2.2x2.0kpc), and an rms 
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Fig. 2.— CO(J=2^1) emission line map of BRI 1335-0417 at 
~1.0kpc resolution, imaged using B array data only and robust 
weighting. Contours are shown at (—3, -2, 2, 3, 4, 5, 6)x<7 
(lcr = 58 fijy beam -1 ). The beam size (0.16 " xO.14 ") is shown 
in the bottom left corner. The labels indicate the rest-frame peak 
brightness temperatures Jj, in the brightest molecular line emission 
regions. 

of 150 /^Jy beam -1 per channel. These channel maps are 
combined into a three color overlay in Fig. [4] 

3. RESULTS 

3.1. Morphology of the Molecular Gas Reservoir 

In Figure[l] the integrated CO(J=2^1) emission over 
the full measured bandpass (308 kms -1 ) is shown at 
a linear resolution of ~1.4kpc (0.2"). The emission 
is clearly spatially resolved over many beams, and ex- 
tended out to a scale of 5kpc (~0.7", deconvolved for 
the beam size). Two distinct peaks of the emission are 
identified, both of which are extended, and embedded in 
a continuous, morphologically complex molecular struc- 
ture. The northern component has a peak flux den- 
sity of 216 ± 50//Jy beam -1 , and the southern com- 
ponent has a peak flux density of 388 ± 50/iJy beam -1 
(peaks labeled 'N' and 'S' in Fig. HJ. In Figure U a 
zoomed-in version of this map at 1.0 kpc (~0.15") lin- 
ear resolution is shown. At this resolution, some of 
the more diffuse structure is resolved out. Clearly, the 
southern component breaks up into multiple line emis- 
sion regions. Most notably, a resolved subcomponent 
extends from the main southern peak ('S') toward the 
northern peak ('N'). In the higher resolution map, the 
northern component has a peak flux density of 246 ± 
58/zJy beam -1 , and the southern component has a peak 
flux density of 375 ± 58 fj,3y beam -1 . The connect- 
ing component has a peak flux density of 202 ± 58 /iJy 
beam -1 . These three peaks alone encompass a molecular 
gas mass of M(H 2 ) = 4.1 x 1O 1O M , corresponding to 

7 As CO emission in currently known z > 4 quasar host 
axies appears to be thermalize d typ ically up to J=4— >3 



(Riechers et al. 2001; Weifi et al. 2007), constant T b is as 
sumed between CO(J=2^1) and CO(J=1->0). The CO is 
more widespread than in nearby ultra-luminous infrared galax- 
ies (ULIRGs), but shows similar physical properties. We thus 
adopt a low ULIRG CO luminosity to H2 mass con version factor 
of c*=0.8Mq (K kms" 1 pc 2 ) -1 (Downes & Solomon [1998) rather 
than «=4— 5 Mq (K km s _1 pc 2 ) -1 as in nearby spiral galaxies 
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Fig. 3. — Channel maps of the CO(J=2— >1) emission in BRI 1335-0417. The same region is shown as in Figure [TJ One channel 
width is 6.25 MHz, or 44 kms -1 (frequencies increase with channel number and are shown at 42615.050, 42621.300, 42627.550, 42633.800, 
42640.050, 42646.300, and 42652.550 MHz). Contours are shown at (-3, -2, 2, 3, 4, 5)xct (la = 150 fj,Jy beam -1 ). The beam size (tapered 
to 0.32 "xO. 30") is shown in the bottom left corner; the crosses indicate main peaks of the the northern ('N') and southern ('S') CO 
components in Fig. [2] 



about two thirds of the mass of Af(H 2 ) = 6.6 x 10 10 M 
derived from the integrated emission seen in the natu- 
rally weighted dataset (7co(2-i)=0.43±0.02 Jy kms" 1 , 
in agreement with previous results by Carilli et al. l2002f ). 
Accounting for the flux in the linewings that is not cov- 
ered by our observations, the full reservoir has an esti- 
mated mass of M tot (H 2 ) = 9.2 x 10 10 M . 

At z=4.4074, the brightest peaks of the northern and 
southern components in the highest resolution map cor- 
responds to beam-averaged, rest-frame brightness tem- 
peratures of T b =61 ± 9K and 40 ± 9K (see Figure [2]). 
Averaging over the area of the full CO( J=2— >1) reservoir 
gives Xb=10.0±0.5K, which implies that part of the dif- 
fuse molecular structure may be rather cold. Note that 
the minimum excitation temperature of CO(J=2^1) 
is T CX =16.6K, which however is only 1.8 K above the 
cosmic microwave background (CMB) temperature at 
2=4.41. 

3.2. Dynamical Structure of the CO Distribution 

In Figure [3J the CO(J=2— >1) emission is shown in 
seven 44 kms -1 wide velocity channels, smoothed to a 
linear resolution of ~2 kpc. The emission is clearly dy- 
namically resolved and is moving, to first order, from 
north to south between the red and blue velocity chan- 
nels (i.e., with increasing channel number). In addition, 
the center channel (#4) shows a bright, compact peak, 
which corresponds to a narrow velocity component. This 
peak is weighted down in the integrated line maps, and 
corresponds to the third small 3cr peak of the southern 
component in Figure (labeled with its lb of 29 ± 9 K). 
Figure [4] shows a composite color map of the velocity 
channel maps and a position-velocity (p — v) diagram 
(left panel) . The middle panel includes all velocity chan- 
nels, while the central channel (covering the range from 



(e.g., Scoville Sz Sanders fl987l; Solomon & Barrett [1991). Such 
low a are also found for z ~ 2.5 submm galaxies (Tacconi et al. 
I2008T) . 



-22 kms -1 to +22 kms -1 ) is excluded from the right 
panel. The emission is clearly moving from redshifted 
to blueshifted emission in the right panel. However, the 
emission in the (green) central channels does not follow 
the north-south extension, but also show an east-west ex- 
tension, indicating an additional, distinct velocity com- 
ponent in the system. The overall peak of the emission 
where all components overlap (seen as a 'white spot' in 
the figure) corresponds to the southern peak in the in- 
tegrated line map. The velocity structure thus is more 
complex than that of a simple, inclined, rotating disk 
(as also seen in the p — v slice). The central, compact, 
narrow velocity component is likely associated with the 
optical quasar. However, the relative astrometry of the 
radio and optical observations (Storrie-Lombardi et al. 
Il996h is unfortunately not accurate enough at present to 
substantiate this conclusion. 

4. DISCUSSION 

We present maps of the molecular gas distribution in 
a quasar host galaxy at z=4.4 at 1.0 kpc spatial resolu- 
tion. The VLA data show that the molecular gas reser- 
voir in this galaxy is not only distributed over a scale 
of ^5 kpc, but also structured in velocity space. This 
is the first time that the molecular gas in a quasar host 
galaxy at such a high redshift (or, indeed, at any red- 
shift greater than ~0) has ever been resolved both spa- 
tially and dynamically over multiple beams (except for 
the le nsed system PSS J2322+1944 at z=4.12; Riechers 
et al. I2008D . The emission in BRI 1335-0417 is resolved 
into at least three distinct components harboring at least 
1-2 x 10 10 Mq of molecular gas each, and at least an- 
other 2 x 10 10 M Q is found in the more diffuse molecu- 
lar medium in between these concentrations. Each of the 
subcomponents hosts a few times the molecular gas mass 
of ne arby ULIRGs such as Arp 220 (Downes & Solomon 
1998). The brightest peak corresponds to a lkpc diam- 
eter region with a gas surface density of E^ c 9 ak =2.3 x 
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Fig. 4. — RGB composite color map representation and position-velocity (p — v) diagram of the CO(J=2— >1) velocity channels shown 
in Fig. [3] with three colors encoding the velocity range of the emission [red: redshifted (channels 1-2 in Figure [3J, green: central (channels 
3-5), blue: blueshifted (channels 6-7)]. The velocity range covered (relative to the line center) in kms - 1 is indicated in the top left corner, 
and the linear scale is indicated in the bottom left corner of each panel. Left: p — v diagram along a position axis of 189° , connecting the 
two brightest peaks of the integrated line emission. Middle: All velocity channels (1-7). The dashed line indicates the orientation of the 
p — v slice. Right: All except the central channel (channel 4). 

source dominantly originates from extended, kpc-scale 
regions and sh ows a peak Tb of only few times 10 5 K 
(Momjian et al. l2007f ). Such Tb are typical for relativistic 
electrons from supernova remnants and the interstellar 
medium in nuclear starbursts, but by at least 2 orders 
of magnitude lower than found in AGN-powered envi- 
ronments. This finding is also reflected in the fact that 
BRI 1335-0417 follows the radi o-FIR correlation for star- 
forming galaxies (Carilli et al. [1999). The FIR-derived 
SFR is by about an order of magnitude higher than in 
local ULIRGs; however, the ratio of SFR and total gas 
mass is consist ent w ith the scaling relation for ULIRGs 
(Solomon et al. 119971) a nd other high-z galax ies (Solomon 
& Vanden Bout 120051 : Riechers et al. I2006D . Moreover, 



10 10 Mq kpc 2 , which is by a factor of a few higher than 
in in the ce ntral regions of typical ULIRGs (e.g., Wilson 
et al. 2008). The total molecular gas mass and distribu- 
tion are reminiscent of those in the z=4.69 and z=6.42 
quas ars B R 1202-0725 and S PSS J1148+52 51 (O mont 
et al. [I996t Ca rilli et al. [20021 Riechers et al. [20061 ; Wal- 
ter etal. " 12001 . albeit concentrated on almost an order 
of magnitude smaller scales than in BR 1202-0725. The 
brightest peak of the CO(J=2^1) emission traces a re- 
gion of 1 kpc diameter where the gas has an average tem- 
perature of at least 61 K (depending on optical depth and 
beam dilution), which is at the high end of but compa- 
rable to temperatures in the centr al reg ions of nearby 
ULIRGs (e.g., Downes & Solomon [1 998). However, in 
ULIRGs, regions with molecular gas temperatures of 30- 
60 K are much more compact than lkpc. The CO exci- 
tation ladder in BRI 1335-0417 is not constrained well 
at present. However, the excitation so far is consistent 
with being similar to the z=4.69 quasar host galaxy of 
BR 1202-0725, for which models of collisional excitation 
predict a kinetic gas temperature of Tkin=60 K and a me- 
dian gas density of n(H2) = 10 cm~ 3 (Riechers et al. 
[2001) . If indeed T kin ~ TP cak , this would indicate that 
the brightest CO peak is close to being resolved. Also, it 
appears that the gas temperature and density are on the 
high end of but comparable to ULIRGs, while the total 
gas mass is by about an order of magnitude higher. 

Based on the FIR luminosity and SED shap e of the 
source (Lfir^S.IxIO 13 ^; Benford et al. 1999), a star 
formation rate (SFR) of 4650M Q yr _1 and a total dust 
mass of 2xl0 9 M Q can be derived. 8 This does not ac- 
count for possible heating of the dust by the AGN. How- 
ever, if the dust in this source has an extension and 
non-uniform distribution similar to that of the molec- 
ular gas, local heating (rather than heating of a cen- 
tral source) is likely to be responsible for most of the 
dust emission. This is consistent with the finding that 
the (rest-frame) 4.0 cm radio continuum emission in this 



1998) 



Assuming SFR=1.5xlO~ 10 L Fm (M e yr -1 / Lq) (Kennicutt 



in the centers of nearby ULIRGs, £P cak correlates with 

r k 

the SFR (Wilson et al. l2008h . which suggests that an in- 
creased availability of fuel for star formation at a certain 
density leads to an increased SFR. Given its high S^ c 2 ak 
on kpc scales, this also motivates the high derived SFR 
for BRI 1335-0417. Assuming that the gas is converted 
into stars at an efficiency of 5 -10% as in giant mo lecula r 
cloud cores (e.g., Myers et al. 119861: Scoville et al. ll987f) . 
the SFR corresponds to a gas depletion timescale of 2-4 
x 10 s yr. 

If the molecular gas in this system was gravitation- 
ally bound, the CO linewidth (see Sect. 2) and distri- 
bution would predict a dynamical mass of M^ yn =l.O x 
10 11 sin~ 2 iMq, which could account for both the to- 
tal molecular gas mass and the mas s of t he black hole 
of M B h=6 x 10 9 M Q (Shields etal. 120061) . but not for 
a substantial fraction of a ^4xlO 12 M stellar bulge as 
predicted if the l ocal M rh-j,, relation were t o hold (Fer- 
rarese & Merritt [200(1 Gebhardt et al. l2000l ). 9 

However, the overall structure of BRI 1335-0417 looks 
rather disturbed, both spatially and in its velocity struc- 
ture. While the general north-south extension of the 

9 In principle, assuming an extreme inclination toward the line- 
of-sight would significantly increase M<j yn ; however, the implied 
large CO linewidths suggest that such a dynamical structure would 
not be virialized (Riechers et al.[2008). 
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source may be in agreement with rotating structure, this 
is not the case for the central part of the emission, in 
particular the compact peak seen in the central chan- 
nel. This, and the structure seen in the high-resolution 
CO( J=2— s-1) map (Figure[2]), are more reminiscent of the 
disturbed gas reservoirs in major mergers, such as seen 
in the nearby Antennae (NGC 4038/39, e.g., Wilson et al. 
2000). If BRI 1335-0417 was an interacting system, the 
northern component may be in the process of merging 
with the southern component, which likely hosts the lu- 
minous quasar. The connecting component then may be 
the part where the two galaxies overlap and merge. Note 
that this region alone would be massive enough to host 
more than 10 of the largest molecular complexes f ound i n 
the overlap region of the Antennae (Wilson etal. 2000). 

Cosmological simulations predict that the merger rates 
at z—4A are substantially higher than at z=0. Such sce- 
narios imply that during the early epoch of hierarchical 
galaxy formation, some of the most massive galaxi es form 
in major merger events (e.g., Springel et al. l2005h . Such 
major, 'wet' mergers are believed to commonly trigger 
both AGN and starburst activity, and lead to high exci- 
tation of the molecular gas during both the hierarchical 
buildup of the h ost g alaxy and the quasar phase (e.g., 
Narayanan etal. 120081 ). In simulations, such objects of- 
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ten show multiple CO emisson peaks, arising from molec- 
ular gas concentrations that have not yet fully coalesced. 

We conclude that the observed properties of BRI 1335- 
0417 (AGN and extreme starburst activity, high CO exci- 
tation, disturbed morphology over 5 kpc scales) are con- 
nected to the ongoing buildup of the quasar host galaxy. 
Such a signpost of early galaxy assembly then could 
be considered direct observational proof of the scenar- 
ios proposed by cosmological simulations, and enable us 
to directly investigate the connection between quasar ac- 
tivity and high-mass merger events at early cosmic times. 
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